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We construct strangelets, taking into account electrostatic effects, including Debye screening, and 
arbitrary surface tension o of the interface between vacuum and quark matter. We find that there 
is a critical surface tension acrit below which large strangelets are unstable to fragmentation and 
below which quark star surfaces will fragment into a crystalline crust made of charged strangelets 
immersed in an electron gas. We derive a model-independent relationship between acrit and two 
parameters that characterize any quark matter equation of state. For reasonable model equations of 
state, we find ffcrit typically of order a few MeV/fm . If cr ^ crcrit, the size-distribution of strangelets 
in cosmic rays could feature a peak corresponding to the stable strangelets that we construct. 
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I. INTRODUCTION 

Ordinary matter consists of atoms, whose nuclei are 
ultimately composed of up and down quarks. One can 
think of nuclei as droplets of nuclear matter, which is 
observed to be very stable: the most stable nuclei have 
lifetimes longer than the age of the universe. However, 
it has been hypothesized [H Q, Q that nuclear matter 
may actually be mctastablc, and the true ground state of 
matter consists of a combination of up, down, and strange 
quarks known as "strange matter" . Small nuggets of such 
matter are called "strangelets" . If this "strange matter 
hypothesis" is true, then many if not all compact stars 
are not neutron stars but "strange stars": large balls 
of strange matter (for a review see Ref. Q). Moreover, 
strangelets may be produced in compact star collisions, 
and hence contribute to the cosmic ray background, or 
in heavy-ion collisions, where they could be observed in 
terrestrial experiments. 

Until recently, it has been assumed that the boundary 
between strange matter and the vacuum is a simple sur- 
face, with a layer of positively charged quark matter be- 
neath it and with electrons floating above it, sustained by 
an electric field which could also support a small normal 
nuclear crust in suspension above the quark matter Q, 
as long as the strange star is not too hot p. 

This was recently questioned by Jaikumar, Reddy, and 
Steiner Q , who showed that if Debye screening and sur- 
face tension were neglected then the surface must actu- 
ally fragment into a charge-separated mixture, involving 
positively-charged strangelets immersed in a negatively 
charged sea of electrons, presumably forming a crystalline 
solid crust. At the surface of this strangelet-|-electron 
crust, there would be no electric field. This gives a radi- 
cally different picture of the strange star surface, making 
it much more similar to that of an ordinary neutron star, 
and casts doubt on all of the previous work on the phe- 
nomenology of strange star surfaces. 

In this paper we address the stability of electrically 
neutral bulk quark matter with respect to fragmenta- 



tion into a charge-separated strangelet-|-electron crust, 
and the almost precisely equivalent question of whether 
large strangelets are stable with respect to fission into 
smaller strangelets. We shall include energy costs due to 
Debye screening, neglected in Ref. Q, and surface ten- 
sion. Our treatment of Debye screening follows that of 
Hciselberg Q except that we include the energy gained 
by charge separation which was neglected in Ref. |31 and 
whose importance was discovered in Ref. |7|. 

We identify the three microscopic parameters that play 
a crucial role in the stability analysis: (1) the charge 
density of quark matter at zero electric charge chemical 
potential denoted as ng; (2) the electric charge suscepti- 
bility of quark matter denoted as XQ'i and (3) the surface 
tension between quark matter and vacuum, denoted as a. 
In Sections ITU and Hill we derive all of our results in terms 
of these parameters. In Section^ we give a parametric 
estimate of the critical surface tension below which neu- 
tral quark matter and large strangelets are unstable to 
fragmentation, and in Section IIIII we construct explicit 
strangelet profiles and thus determine the dimensionless 
quantities left unspecified in Section ITU 

The quantities uq and XQ can be evaluated for any 
proposed quark matter equation of state, and in Section 
II VI we shall do so in bag models for unpaired quark mat- 
ter and quark matter in the 2SC color superconducting 
state 13, [13, llil [l3l- We find that even if the strange 
quark is so heavy that no BCS pairing occurs in bulk 
quark matter, if such bulk matter is unstable to fission 
then the strangelets that result arc likely to be in the 
2SC phase. This occurs because charged 2SC matter 
may have higher pressure than charged unpaired quark 
matter, even when neutral 2SC matter has lower pressure 
than neutral unpaired quark matter. 

The color-flavor locked (CFL) phase [13,113 is the fa- 
vored phase of neutral quark matter if the density is high 
enough, the pairing interactions arc strong enough, or the 
strange quark is light enough. The CFL phase is an insu- 
lator, so both nq and XQ vanish, because there is pairing 
between all three flavors, and no electrons are needed in 



order to maintain neutrality. This makes CFL matter 
stable against fission even if a is arbitrarily small. Our 
goal in this paper is to understand what happens if the 
quark matter at the surface of a strange star is not in 
the CFL phase, although this form of quark matter may 
occur deeper within. 

We close in Section with comments on phenomeno- 
logical implications for the strange star surface and ex- 
perimental strangelet searches. 



II. STABILITY OF STRANGELETS 

The stability of strangelets is determined by the equa- 
tion of state (EoS) for quark matter (QM), which will 
be discussed below, and the EoS of the vacuum with 
electrons, which is well known. The equation of state 
expresses the pressure p(/i, /Ue) as a function of the chem- 
ical potentials for the two conserved quantities, quark 
number (/i) and negative electric charge (/ie)- We work 
at zero temperature throughout. The quark density n 
and the electric charge density q (in units of the positron 
charge e) are 



n = 



dp 



q 



p 



dp 



(1) 



The chemical potential fi^ is the energy cost of adding an 
electron, so it is the positron charge e multiplied by the 
electrostatic potential cj). Hence, /ie9 = (f'P- The energy 
densities follow from the usual thermodynamic relation 



/in — fi^q ^ p 



(2) 



The second term in (|2} and the final term in ^ arc 
negative because fie is the chemical potential for negative 
electric charge. As written, e and p are understood to 
include the electric field energy and pressure. Henceforth, 
however, we shall intepret p as only the kinetic pressure of 
the quark matter and electrons, in which case the electric 
field energy i J (Pr(f>p must be added to / (Pre. This can 
be implemented by writing 



fj,n - 



jMe9 -PQM+c 



(3) 



We make the strange quark matter hypothesis that 
at some fj-cnt, uniform neutral quark matter with some 
nonzero n has the same pressure as the vacuum. (The 
hypothesis further requires that /Xcrit < Mnuciear, where 
Mnucicar ~ 310 McV is the chemical potential of uniform 
neutral nuclear matter at zero pressure, and that the 
chemical potential for two-flavor quark matter with zero 
pressure is greater than /inucicar-) However, if only vol- 
ume contributions to the energy are taken into account, 
then low pressure quark matter is unstable against frag- 
mentation into positively charged strangelets embedded 
in a gas of electrons [3. In this paper we evaluate the 
surface tension required to stabilize neutral quark matter 
at zero pressure against such fragmentation, taking into 
account Debye screening. 



Debye screening causes the positive charge density in 
large strangelets to migrate towards the surface, resulting 
in a charged skin whose thickness is of order the Debye 
length Xu and a neutral interior. So, for large strangelets 
{R 3> \d) the volume energy benefit of charge separa- 
tion is reduced to a surface energy. Typically, Xd is of 
order 5 fm in quark matter. Whereas Debye screening pe- 
nalizes large strangelets, surface tension penalizes small 
strangelets most, as their surface to volume ratio is the 
greatest. Both effects make fragmentation into positively 
charged strangelets less favorable. We can expect their 
combined effect to be least for droplets with R ~ A^, 
meaning that if bulk quark matter proves unstable, it 
will fragment into strangelets of this size. 

Jaikumar et al. proposed that the outer layers of a 
strange star could be constructed from positively charged 
strangelets with some radius R embedded in a gas of 
electrons, forming a Wigner-Seitz lattice |3- The local 
pressure p will vary considerably within a single Wigner- 
Seitz cell, but the external pressure Poxt on a Wigner- 
Seitz cell from its neighbours will be zero at the surface of 
the star, and increase smoothly with depth. To evaluate 
the stability of such a crust, we must evaluate the Gibbs 
free energy per quark 



9iR) = 



E{R)+p,^tV{R) 
N{R) 



(4) 



where E, V and N are the energy, volume and quark 
number of a Wigner-Seitz cell containing a strangelet of 
radius R. We shall calculate the energy E{R) by solv- 
ing the Poisson equation upon making a Thomas-Fermi 
approximation to obtain the distribution of charge, inte- 
grating e as defined in Q. and then adding the surface 
tension energy AnR^a. We must then compare g{R) with 
the Gibbs free energy per quark ggM of uniform neutral 
quark matter at the same pressure Pext- We shall work 
throughout at pcxt = 0, meaning that we analyze the 
stability of the outer surface of the mixed phase crust 
proposed in Ref. [3 and at the same time that of an iso- 
lated strangelet with radius R. We shall be comparing 
them to neutral bulk quark matter with /i = /^crit, which 
has zero pressure and ggM = /icrit- So, at poxt = 0, the 
free energy per quark of a strangelet of radius R, relative 
to infinite neutral quark matter, is 



Ag{R) = giR) - gQM 



E{R) 

N{R) 



/^crit 



(5) 



For any a, Ag{R) — > for R ^ 00. For large enough 
values of a, Ag(R) will be dominated by its surface en- 
ergy contribution 3a/{nR), making Ag(i?) > for all R. 
This means that for large enough a neutral bulk matter 
is stable with respect to fragmentation. For large enough 
a, furthermore, Ag{R) decreases monotonically with in- 
creasing R, meaning that isolated strangelets of any size 
are stable, but can lower their energy by fusing with other 
strangelets should they encounter them. At the other ex- 
treme, with a — wc know from the work of Ref. [3 that 



Ag{R) < for small enough R. We therefore expect that 
for small enough values of a we shall find a range of R in 
the vicinity oi R ^ X^ for which Ag(i?) < 0. This means 
that, for small enough a, neutral bulk quark matter and 
large strangelets are unstable to fragmentation, with the 
stable strangelets being those having the size R = R^ 
that minimizes /S.g{R). 

The equation of state of quark matter at phenomcno- 
logically interesting densities and T ^ fi cannot cur- 
rently be calculated by lattice gauge theory or by other 
methods, so we can either obtain an approximate EoS 
from some model, or make a general parameterization. 
We will find that /Xe ^ /x for all strangelets in all models 
that we consider, so a general parameterization of the 
EoS can be obtained by expanding in powers of /^e/A* |j9; 

PQM ^Po{n,ms) -nQ{fi,ms)fie + ^XQ{^J','rns)^il + . . . 

(6) 
This second-order expansion, which neglects the electron 
pressure p^ ~ /x^, can be used for any model EoS or for 
that predicted by QCD. We will sec below that it is an 
excellent approximation for the analysis of strangelets. 
It reduces the EoS-dependence to specifying the three 
functions po, nq and XQ- Moreover, only nq and XQ 
occur in the Poisson equation. 

Neutral bulk quark matter has g == 0, so its electron 
chemical potential ^ is 



/^e 



utral _ ^ 
~ XQ 



(7) 



The quark chemical potential /icrit at which it has zero 
pressure is determined by solving 



PO (Merit, "Is) 



nniflcrit^ms) 



(8) 



2XQ(/^crit,ms) 

The Debyc screening length in quark matter is controlled 
by XQ, and is given by 

1 



Xd = 



y^AnaxQ 



(9) 



with a = 1/137. 

We close this section with a parametric estimate of the 
critical surface tension below which strangelets with R ^ 
Xd have Ag{R) < 0, making neutral bulk quark matter 
(and larger strangelets) unstable to fragmentation. In 
section IlIII we will quantitatively determine /Xe(7') and 
the radius i?* of the most stable strangelet. For now, 
we write x« = R.,/Xd and take /Xg in the strangelet to 
be a constant with the value Juq/xq- Presuming the 
dimensionless quantities .t* and / to be of order unity, 
we use Eqs. (|TJ, ®, ©, ®, O, ® and I® to obtain 



A.g(i?,) 



3(7 



'Q 



(1-/) 



ni?* 2riXQ 

which is negative if cr < cTcrit where 

_ (1 - f)R.nl _ (1 - f)x,nl 



6XQ 



/ 3/2 

/Traxg 



12, 



(10) 



(11) 



In the next section we shall construct strangelets by solv- 
ing the Poisson equation and show that this parametric 
estimate is valid, with x^ w 1.61 and / w 0.49. 



III. CONSTRUCTING STRANGELETS 

We assume that the lowest energy state of a strangelet 
of radius R is spherically symmetric. The quark chemical 
potential ^ is independent of r because the only net force 
on a given part of the strangelet is the electrostatic force. 
(We assume our strangelets are small enough that their 
self- gravity is unimportant, and that color is screened so 
strong interactions do not occur across distances greater 
than about 1 fm.) The value of pi inside the droplet 
is a little higher than /icrit because the surface tension 
compresses the droplet slightly. To determine the value of 
/i, we require the pressure discontinuity across the surface 
of the strangelet to be balanced by the surface tension: 



lim (piR 



Sr) - p{R + Sr)) 



2fT 



(12) 



The larger a becomes, the more compressed the 
strangelet is, and the higher the value of /i inside the 
strangelet. 

To calculate fieir) — e(j){r) in the Thomas- Fermi ap- 
proximation we solve the Poisson equation, which in 
Heaviside-Lorentz units with h = c = 1 takes the form 



V20(r) 
i.e. V^fieir) 



-p{r) , 
-A-Kaq{r) 



(13) 



subject to the boundary conditions 
lim r(j){r) ~ 



Zr, 



-f(r = 0) = 
dr 

^(i?-)-^(i?-) = o 

dr dr 



(14) 



The second and third boundary conditions follow from 
the fact that there are no delta-functions of localized 
charge, so (/)(r) and d4>/dr are continuous everywhere. 
The first boundary condition states that the net charge 
of the strangelet, including any electrons inside or around 
it, is Zooe. The Fermi wavelength of the quarks is 2-k / pi 
which is around 4 fm, so it is reasonable to use the 
Thomas- Fermi approximation to describe the charge dis- 
tribution due to quarks inside a strangelet with diameter 
of order 3 Ad ~ 15 fm, which is the approximate size of 
the stable strangelets that we will find. 

Strangelets with size R ~ Ad are too small to have 
electrons localized within them, meaning that the charge 
of the strangelet itself is given by that of the quark mat- 
ter, 



Z = 



d rq{r) = / d r {nq - XQfJ'ei.r)) 



(15) 



To analyze an isolated strangelet, not surrounded by an 
atom-sized cloud of electrons, we must find solutions 
to the Poisson equation with Zoo = Z. To analyze 
a strangelet located within the Wigner-Seitz cell of a 
strange star crust, however, we must include a cloud of 
electrons around the strangelet such that the strangelet 
and electrons together satisfy the Poisson equation with 
^cx) = 0. (Each Wigner-Seitz cell contains one strangelet 
and Z electrons, ensuring that the crust is neutral on 
macroscopic length scales.) In our numerical work, we 
have included the electrons outside the strangelet, tak- 
ing qir) = — /ig/37r^ for r > R, and found solutions to 
the Poisson equation with Zoo = 0. As in atomic physics, 
the Thomas-Fermi description of the cloud of electrons 
around a strangelet with charge Z becomes more accu- 
rate with increasing Z. Fortunately, it turns out that the 
electrons do not play an important role in the stability 
of the strangelet, so we can ignore them, setting q{r) = 
for r > R and finding solutions to the Poisson equation 
with the boundary condition at infinity given simply by 
Zoo = Z. We have done all our analytic calculations with 
this boundary condition, and have confirmed numerically 
that adding electrons to the system in order to find so- 
lutions with Zoo = makes a negligible difference to the 
energy E and to fJ.e(r) inside the strangelet. This sim- 
plification arises because, as mentioned in Section II, the 
electron contribution to the EoS is subleading in fie/ fJ^- 

We obtain E{R) and N{R), the quantities which ac- 
cording to (0) we need in order to determine the stability 
of a strangelet with some specified radius R, as follows. 
For a given fi we solve (|13|l subject to the boundary con- 
ditions (|14|l . We repeat this for different values of /i until 
we find the one that obeys l|12|l . We then have fi and 
fieii"), so from ^ and (^ we can obtain the energy den- 
sity e(r) and quark number density n{r). We integrate 
these, and add in the surface energy, to obtain the total 
energy E and quark number N and hence the energy per 
quark E/N. 

The solution to the Poisson equation satisfying the 
boundary conditions (|I4|I with Zoo = Z given by H15|) 
is H 



r<R: fie{r) = 7^ ( 1 

r>R: Heir) 
where 



riQ I Xp sinh(r/Aji) 
XQ \ r cosh{R/\D) 

Za 



Z = A-kuqXj^ {x — tanhx) , 



(16) 



(17) 



where we have defined x = R/Xd- We see that /ie(r) < 
^neutral „ fiQ / XQ throughout the strangelet, with fJ,e{r) 
closer to zero for smaller strangelets and closer to ^'^™^''^^ 
for larger ones. 

To this point we have not determined jj,. (And, recall 
that both uq and XQ depend on /z.) Because we have 
no electrons and hence p = outside the strangelet, the 
condition (|12|l requires that just inside the surface of the 




FIG. 1: A strangelet profile with radius R = I.6X0, and 
a = 0.1325 nQXo/XQ- (This corresponds to the barely-stable 
strangelet at the critical surface tension — see Section lllll l 
The horizontal axis is 2; = r/Xo- The charge density q is 
plotted in units of nq, /j.^ is plotted in units of uq/xq, the 
electric field is plotted in units of nqj ^fxQ a-nd the pressure 
Pqm is plotted in units of n^/xQ- Because of the electric 
field, the charge density in the strangelet is pushed towards 
the outer edge. The pressure gradient within the strangelet 
is balanced by the electric force on the charged matter. 



strangelet the pressure is given by 2cr/R. Recalling that 
n =: dp/dfi and that p{p.c-cit) = and using Eqs. © and 
(I16|l we see that this requires that the quark chemical po- 
tential inside the strangelet (which is constant through- 
out its interior) is given by 



M — /^crit 



'Q 



tanh^ : 



2"-XQ 



2cr 



^ ■ (18) 



The negative term present even when cr = arises be- 
cause p.e{R) < ^|;j™'™i, making p(/^crit, Me (-R)) > 0. To 
achieve p = 0, as required just inside the edge of the 
droplet if CT = 0, /i must be less than /^crit- The a- 
dependent positive term then enforces the pressure dis- 
continuity 1)12(1 required at nonzero a. (The derivation 
requires that each of the corrections to fi are separately 
much smaller than //crit, which is well satisfied in all re- 
sults we show.) 

We can now construct the profiles of fJ,e{r), q{r), the 
electric field and p{r) for strangelets with radius R. An 
example is shown in Fig. ^ We show the profiles in 
terms of hq and XQ- Expressed this way, they are model- 
independent. Given a model equation of state, hq and 
Xq must be evaluated at the /i given in (|18|l . which itself 
depends on ng and XQ- This means that the explicit 
determination of /x in a model must be done numerically, 
although in practice we find that the critical surface ten- 
sion is very small, so for a « iTcrit we can use nq {p,crit) 
and XQifJ-ciit) without making significant errors. We have 
compared the profiles obtained analytically in terms of 



uq and XQ ^s we have described with those obtained 
numericaUy, with ng and XQ specified by the bag model 
expressions for unpaired quark matter given in Section llVl 
and with electrons included in the solution of the Poisson 
equation. The agreement between analytical and numer- 
ical profiles is excellent. 

We now evaluate the Gibbs free energy per quark of 
a strangelet of radius R relative to that of bulk neutral 
quark matter, namely ^g{R) of (|3J). Integrating the en- 
ergy density from Q and adding the surface contribu- 
tion, E{R) is given by 

/■^ 
E{R) = / d^r [-p{r) + ^n - ig(r)^e(?')] + '^nR'^a 
Jo 

d^r [-pa + fin+ ^72Qfie{r)] + i-KR^a , 

(19) 

where wc have used ^ and © to replace p{r) and q{r). 
We see now that the only property of the profile fJ-eif) 
that we need in order to evaluate E{R) is the volume 
average 



-ft Jo XQ XQ 



tanhx 



(20) 

where / is the parameter we introduced in making the 
estimates (|10|l and (|ll|l . and which we have now evalu- 
ated using the profile fJ,e{r) in H16|l . Thus, by solving the 
Poisson equation we have learned that 

X — tanh X , , 

l-/ = 3 . . (21) 



From (UD) and ^^ we find 



N{R) - 
and hence, using 



Po 
n 



"^nxQ 



and (O, 



nAg{R) 



2XQ 



Po +n- (m- Merit)- 



3cr 
nR 



)(1-/) 



(22) 



3cr 

2xQ ~R ■ 

(23) 
Notice, however, that the first term is zero when evalu- 
ated at /i = /icrit and is given by n(/Xcrit — m) to lowest 
order in (fi — /icrit)- Neglecting fractional errors of order 
(/i — /icrit)//icrit in Ag (which Fig.|21shows are negligible) 
we obtain 



nAg{R) 



3nj 



tanh : 



2XQ 



3ct 

'r 



(24) 



which is (|1U|I with (1 — /) now known but R* still to 
be determined. (By comparing to solutions obtained nu- 
merically, we have confirmed that at large values of a/R 
the fractional error introduced in Ag is ~ a / (nRficiit) ^ 
as in (|18|) .) It is convenient to write (|24|l in terms of a 




x=R/Xd 



FIG. 2: The Gibbs free energy (i.e. energy per quark) relative 
to neutral uniform quark matter with p = for strangelets 
of various sizes and surface tensions. We plot Ag{x), which 
is Ag in units of nQ/{nxQ), as a function of the radius 
X = R/\d, for various values of the surface tension (a is 
a in units of Wq Ad/xq)- The solid lines (red online) are ob- 
tained from (I24II : the dashed lines (blue online) are obtained 
from a numerical solution of the Poisson equation, includ- 
ing the electrons and not making any approximations in the 
evaluation of E{R). The numerical results confirm the va- 
lidity of the approximations made in deriving I|24I25|I . (The 
numerical analysis is for a bag model of unpaired quark mat- 
ter with /icrit = 305 MeV and iris = 200 MeV. The nq and 
XQ for this model are given in Section IIVI ) At the criti- 
cal surface tension CTcrit ~ 0.1325 one can just barely con- 
struct a strangelet that is favored over uniform quark matter. 
The radius of this critical strangelet is R, = 1.606 Ao. For 
CT < CTcrit, large strangelets and bulk neutral quark matter 
are unstable to fragmentation. For a > CTno-barrior ~ 0.1699, 
there is no fusion barrier: strangelets of any size can lower 
their free energy by fusing should they encounter each other. 
For CTcrit < CT < CTno-barrior, there are "metastable" droplets, 
whose free energy is greater than that of bulk quark matter 
but which must overcome a barrier in order to fuse. 



dimensionless function Ag{x) of the dimcnsionless radius 
X = R/Xd, 



Ag{x) 



3 a; — tanh x 3a 
2 ^^3 ^~ 



where 



Ag{R)^^^Ag{x), a^ ^^"^ 



nXQ 



Xon^Q 



(25) 



(26) 



In Fig.|2]we plot Ag versus R in dimcnsionless units for 
several values of cr; this is equivalent to plotting Ag{x) for 
several values of ct. The function Ag{x) has a stationary 
zero at x* = 1.606 for CTcrit = 0.1325: this corresponds 
to the marginally stable strangelet at the critical surface 
tension. For W < o'critj the minimum of Ag occurs at 
smaller a;, and the value of Ag at the minimum is neg- 
ative, corresponding to a smaller stable strangelet. For 



o" < CTcrit, neutral bulk quark matter and large strangelets 
can lower their free energy by fissioning, although for 
large enough strangelets there is always an energy barrier 
to fission. In the marginally stable droplet with ct = CTcrit 
and X = x^, H18() simplifies to ^ = ^crit, with the lat- 
ter two terms in this equation cancelling. This means 
that the quark chemical potential in such a droplet is the 
same as that in bulk neutral quark matter. In a droplet 
with size x = x*, the volume average of j2e is reduced 
from its value in neutral bulk quark matter by a factor 
/ = 0.5051. (Note that we plotted Fig. [2 for a = CTcrit 
and X = x*.) The parameters in the estimate (|ll|l are 
now fully determined. 

We see in Fig.|21that at ffno-banicr == 0.1699, the min- 
imum in A.g{x) becomes a stationary inflection point 
at a; = 2.772. For ^crit < ct < CTno-barrior, then, 
there exist metastable strangelets with sizes ranging from 
R = 1.606 Xd to R = 2.772 A^- These are stable against 
fragmentation, and if two of them encounter each other 
there is an energy barrier to their fusion. However, they 
do have higher energy per quark than neutral bulk quark 
matter. For a > ffno-barricr there is no local minimum, 
and all strangelets can lower their free energy by fusing 
with other strangelets should they encounter them. 

Converting back to dimensionful quantities, we find 



(27) 





nj^Xn 


= 0.1325- 


"Q 


CTcrit 
'^no— barrier 


XQ 
- 0.1699 ^ . 


/i^xf 



XQ 



If the strange matter hypothesis holds (e.g. if /icrit < 
310 MeV) and if a in QCD takes on any value less 
than (Tno-barricr, thcrc wiU bc a favored size for the 
strangelets that could be found in cosmic rays. If 
CT < cTcrit, if any strangelets were found they would all 
have sizes peaked around a single value. If a is some- 
what larger, in the regime where there are metastable 
strangelets, the distribution of strangelet size in the uni- 
verse would include a peak at the size corresponding to 
the metastable strangelets, and a continuous distribution 
of larger strangelets, big enough that they sit beyond the 
local maximum in the Ag{R) curve, where Ag{R) is a de- 
creasing function of R. There would be no energy barrier 
for these larger strangelets to fuse with one another. 



IV. EVALUATING THE STABILITY OF 

STRANGELETS MADE FROM UNPAIRED AND 

2SC STRANGE QUARK MATTER 

To this point, we have not needed to specify a model 
for the quark matter equation of state. Given the values 
of uq and XQ in any such model (or, ultimately, in the 
equation of state of QCD) our results can be used to 
evaluate the critical surface tension dcrit . If cr < acrit in 
QCD, strange stars will have a solid crust formed from 



positive strangelets immersed in an electron gas. And, 
large strangelets will be unstable to fission. 

In this section, we provide two bag model examples 
to illustrate how our model-independent results translate 
into estimates of cTcrit in MeV/fm^. In bag models one 
writes pqm(M7 A'e) = Pqui^iksip-, He)-B where B is the bag 
constant, andpquarks is obtained by making assumptions, 
such as treating the quarks in the bag as noninteracting 
("unpaired quark matter"), or as being in a color su- 
perconducting phase, or perhaps as having some weak 
residual QCD interactions, although we shall not include 
the last effect here. We begin with a bag model in which 
the quarks inside each bag are noninteracting, and then 
discuss the consequences of BCS pairing, which would 
lead to color superconductivity in infinite quark matter. 
We shall set the up and down quark masses to zero, and 
treat the strange quark mass tti^ as a parameter. One ex- 
tension of our work would be to evaluate ng and XQ in 
models in which the (stran ge) qua rk mass(es) are solved 
for self-consistently |14L llSl Il6l Il7| , rather than taken as 
a parameter. 



A. Unpaired quark matter 

A derivation of the pressure in the bag model for un- 
paired, noninteracting, quark matter can be found, for 
example, in Ref . |l8l | . Expanding in powers of /.ig accord- 
ing to © yields 



XQ{f^,ms) 



M 



-B + ^ + — {2^l^ - 5/im3) y/fi^ - ml 
3m4 

37r 



M + Vm^ - -ml 



I 2 2\3/2 



27r2 

2 



h ^^^' + ^^^" '"') ^ 1^ 



(28) 
(29) 
(30) 



In the second expressions for nq and XQ^ ^^ have fur- 
ther assumed that ttIs <^ fJ-. The resulting simplified 
expressions yield results that are familiar in the litera- 
ture, for example the fact that ^^^^^'^'^^ = hq/xq gives 
the well-known lowest-order result ^^«"*'''^i = rn^/^Afj,) 
for unpaired quark matter. We shall not actually use 
these approximate expressions in the following, because 
with n ^ 300 MeV as in a strangelet or at the surface of a 
strange star we cannot necessarily assume that nis <C /i. 

Instead of fixing the bag constant B and varying other 
parameters, we shall fix the more physical quantity /icrit- 
At each value of /icrit and rris we use ^ to fix _B. 

In Fig. 13 we show how dcrit varies with rris in the bag 
model for unpaired quark matter for /icrit = 305 MeV. 
(The curves for the other allowed values of /Xcrit, which 
vary from 283 MeV to 310 MeV, are indistinguishable 
from this one except that, as we discuss below, the nis 
at which they terminate is /icrit-dependent.) We see that 
the maximum value of the critical surface tension is less 




250 



m^ (MeV) 



FIG. 3: The critical surface tension for strangelets made of 
unpaired quark matter, as a function of the strange quark 
mass, at /Xcrit = 305 MeV. Other allowed values of /icrit give 
almost indistinguishable curves. The solid curve is given by 
1271 1 whereas the dashed curve was obtained by solving the 
Poisson equation numerically. 



density n = 0.9093 fm-^. If cr = cr„it = 1.377 McV, the 
critical strangelets with radii I.GOGAd = 7.782 fm have 
baryon number A w 598 and charge Z w 69, consistent 
with the result Z « 0.1^ from the literature [l3. If A 
were significantly smaller, as for example would be the 
case for a significantly below Ucrit , it would be important 
to include curvature energy and shell effects in the cal- 
culation oi E/N. The analysis of Ref. 19| indicates that 
these can reasonably be neglected for strangelets with 
A's greater than a few hundred. 

Finally, let us discuss the effect of interactions. Al- 
ready in the earliest bag model analyses of strange stars 
and strangelets Ql, the perturbative QCD interactions 
between quarks inside a bag were taken into account. To 
zeroth order in ms, these just introduce modifications to 
the relationship between ^ and n which are small if as is 
small. We leave the inclusion of such perturbative effects 
to future work, and turn to the possibility of color super- 
conducting quark matter, which can introduce changes 
to the values of hq and XQ that are qualitative, in the 
sense that they do not become arbitrarily small when 
some parameter of the EoS is varied. 



than 2.7 MeV/fm^. The solid curve in the figure is given 
by H27|) with uq and XQ now specified for the unpaired 
quark matter bag model by (|29|l and (|30|l . The dashed 
curve was obtained by solving the Poisson equation nu- 
merically, including the electrons, and then evaluating 
Ag without making the approximations that went into 
deriving if^ . 

The curves in Fig. O end at m^ = 240 MeV because 
beyond this value nuclear matter becomes unstable rela- 
tive to two-flavor quark matter. Unlike the conversion to 
three-flavor quark matter, which must surmount a barrier 
of many weak interactions, this would allow rapid con- 
version of nuclei to quark matter. For ficHt ~ 305 MeV, 
requiring that the pressure of two-flavor quark matter 
at ^ = 310 MeV be negative requires nis ^ 240 MeV. 
Note that we can consider values of m^ all the way up 
to n, at the expense of tuning ^crit closer and closer to 
Mnuciear = 310 McV as wc take rris — > fi in such a way 
as to keep nuclear matter (just) stable with respect to 
two-flavor quark matter. In this fine-tuned limit, we 
can get strange quark matter with an arbitrarily small 
strange quark density, with uq and XQ arbitrarily close 
to their two-flavor values /i'^/(37r^) and 5/i^/(37r^) re- 
spectively. From (|27|1 we see that in this limit, CTcrit -^ 
0.01185 nl^^^^jVUn^a = 5.502 McV/fm^. This is the 
absolute limit to how large CTcrit can be pushed in the bag 
model for unpaired quark matter. Fig. O shows the limit 
if one requires fJ-ciit to be 5 MeV below /i„ucicar- 

To give a sense of the scales that characterize the crit- 
ical strangelets, let us take noninteracting quark matter 
with /icrit = 305 MeV and m,, = 200 MeV as an exam- 
ple. We find UQ = 0.07104 fm^^ XQ = 0.4644 fm^^ 
and X]j = 4.845 fm. Bulk neutral quark matter has 
^neutral _ Uq / XQ — 30.1 McV and has quark number 



B. Color superconducting quark matter 

The strong interaction between quarks is attractive in 
the color-antisymmetric channel, and this leads to BCS 
pai ring and color superconductivity in cold quark mat- 
ter |l2l| . The critical temperature for color superconduc- 
tivity is expected to be in the range of tens of MeV, 
which is far above the temperature of strangelets or neu- 
tron stars. We therefore work at T = throughout. 

Unlike perturbative interactions, color superconduc- 
tivity can have dramatic qualitative consequences for 
the properties of quark matter, because those quarks 
which undergo BCS pairing have their number densities 
"locked" into being equal, as the pairing energy gained 
in so doing overcompensates for the free energy cost of 
maintaining number densities that would not minimize 
the free-particle free energy. In this section we explore 
the consequences of color superconductivity for the sta- 
bility of strangelets. 

To sketch the context for our discussion, let us describe 
the neutral color-superconducting phases that may oc- 
cur in quark matter with ^ « 310 MeV, depending on 
the strength of the attractive quark-quark coupling that 
leads to BCS pairing. We expect the color-flavor locked 
(CFL) phase at the strongest values of the coupling, 
with some other, probably non-isotropic, phase occurring 
at intermediate coupling, and unpaired quark matter at 
lower couplings still. As we will see, the non-isotropic 
phases are expected to yield stable strangelets similar to 
those arising from unpaired quark matter, whereas the 
analysis of the stability of CFL quark matter strangelets 
is qualitatively different. We will see that the two-flavor- 
paired "2SC" phase is not expected to be the favored 
state of neutral quark matter at ^ w 310 MeV, but it may 
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well play a role in strangelets, which consist of charged 
quark matter. 

The CFL phase [ij, UM i^ ^^^ ground state of three- 
flavor quark matter when the attractive quark-quark cou- 
pling is sufficiently stron g th at its gap parameter sat- 
isfies AcFL > mll{2ij) ill llll. In the CFL phase, 
quarks of all three colors and all three flavors undergo 
BCS pairing [13, and the resultant locking of the Fermi 
surfaces makes it an electromagnetic insulator, neutral in 
the absence of electrons [l^ |23| , with nq = XQ = and 
^neutral _ nq / XQ = 0. If bulk quark matter at ^crit is in 
the CFL phase, strange stars will be neutral with /ie = 
in their interiors |23l | , meaning that there is no possibility 
of charge separation and fragmentation, no possibility of 
a solid crust, and no reason for large strangelets to fission 
into smaller ones. This can be seen from our analysis by 
noting that in the CFL phase iJcrit = 0. 

At intermediate values of the quark-quark coupling, we 
know that a different color superconducting phase must 
occur. For m1/{5.2fi) < Acfl < 'm1/{2fj), model analy- 
ses that are restricted to isotropic phases predict a gap- 
less CFL (gCFL) phase [H HH . Depending on the value 
of Acflm/'^s; the values of nq, XQ and /Uneutrai in gCFL 
quark matter can fall anywhere between those character- 
izing the CFL and unpaired phases, and in particular can 
be arbitrarily small but nonzero. This would make the 
investigation of strangelets in some gCFL-like phase very 
interesting, as the Debye length could be arbitrarily large. 
However, the gCFL phase per se is unstable to the forma- 
tion of current-carrying condensates |2J, |2^ |23, [23, |2l| 
and so it cannot be the ground state. The nature of 
the ground state at intermediate coupling is not yet es- 
tablished, but one possibility is a three-flavor crystalline 
color superconducting phase [23, 123, ISfl |31| , with a non- 
trivial crystal structure like that favored in the simpler 
two- flavor case |32, |33 ■ Such phases do not involve the 
locking of Fermi surfaces: this is in large part why they 
are well-motivated candidates to replace the gapless CFL 
phase. This means that, unlike in the CFL or gCFL 
phases, their values of nq, XQ and hence nq/xQ are 
likelysimilar to those of unpaired three-flavor quark mat- 
ter |23j 13^ , so that strangelets made of such matter will 
have similar stability properties to ones made of unpaired 
quark matter. 

The final possible color superconducting phase we con- 
sider is the "2SC phase" in which BCS pairing occurs only 
between up and down quarks of two colors |3. llOlllll [l3 | . 
Previous studies indicate that we should not expect this 
phase to occur in bulk neutral quark matter near the 
transition to the vacuum (which we are assuming occurs 
at /icrit < 310 MeV). In bag models in which nis is a 
parameter, the neutrality constraint ensures that either 
unpaired, CFL or gCFL quark matter is always favored 
over 2SC at T = [H |2l|. In models that solve for 
nis self-consistently there is sometimes a small 2SC win- 
dow, but only if the interaction and parameters are cho- 
sen such that TTis is much too large for three-flavor quark 
matter to exist with /x = ^crit < 310 MeV [TUTl H a llTJ . 



However, even if the 2SC phase is not favored in 
bulk matter, which must be neutral, it may occur in 
the charged quark matter of the stable (or metastable) 
strangelets that we have constructed. We see in Fig. Q] 
that /^e inside these droplets is significantly less than 
^neutral ^ This mcaus that the unpaired u and d Fermi 
surfaces are closer to each other, and farther from the 
s Fermi surface, than in neutral unpaired bulk matter. 
This should favor the 2SC phase. To investigate this, we 
take results from Ref. [I^l and use them to evaluate the 
parameters of the generic quark matter EoS © , 



nq^{t^,m,) 
Xq'^iP^ms) 



A2 ,,2 
unpaired ^ '-^2SCr^ 



Po 



unpaired 



n 



^2SCA^ 



37r2 



(31) 

(32) 



187r2 
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12/1^ + 6/iv//u2 -ml + A 



2 
2SC 



(33) 



(These expressions have all been derived assuming 
A2SC < M, which implies Assc = 2^/3 Acfl [13 ■) Note 
the change in XQ '■ it is approximately half as large as in 
unpaired quark matter, meaning that the Debye length in 
2SC quark matter is larger than that in unpaired quark 
matter by a factor of V^. This change once again orig- 
inates in the locking of those Fermi surfaces which pair, 
so it is qualitative, in the sense that it is independent 
of the value of A2sc- To see how big an effect this 
could have on the stability of strangelets, imagine giv- 
ing the unpaired phase a larger bag constant than the 
2SC phase, so that bulk neutral quark matter is in the 
2SC phase. The results of section iHll show that the size 
of the critical strangelet would be increased by a factor 
of approximately v2, and tTcrit would be increased by a 
factor of approximately 2-\/2. However, if we keep the 
same bag constant in unpaired and 2SC quark matter 
then the effect is smaller. For example, let us analyze 
the case where m^ = 200 MeV, Assc = 25 MeV and 

^unpaired ^ 3^5 y^^y ^ ^^ g^^ ^2SC ^ 1 . Q 1 2 n^"P"""^ and 

Xrf^ = 0.479 Xq°^^"^° . In this case we cannot naively 



use the results of section IIIII (yielding cTcrit enhanced by 
the factor of about 2>/2 described above) because fi^f^ 

is 305.507 MeV, larger than M™t"^"^'' ' and consequently 
the critical strangelet occurs when the Ag{R) curve con- 
structed for the 2SC phase has a local minimum with 



AgiR) = fJ.'^X - Merit < 0, corresponding to 2SC 



strangelets with the same E/N as bulk neutral unpaired 
quark matter. We then find Ucrit = 2.79 McV/fm^ and 
i?, — 1.16 A|^*-^ ~ 8.1 fm. Comparing this to the re- 
sults for unpaired quark matter strangelets at the same 
/icrit and nis (which gave Ucrit = 1-38 MeV/fm and 
i?» = 7.78 fm), we see that the occurrence of 2SC matter 
has changed CTcrit by a factor of about 2. The increase 
in Xd has been cancelled by a decrease in x* with the 
result that i?* hardly changes. Evaluating A and Z for 
the critically stable 2SC droplet, we find A w 677 and 



Z K 105. 

It should be noted that the 2SC Cooper pairs, whose 
size is ~ I/A2SC, fit within these strangelets, but not 
by much. This means that for a specified interaction 
strength between quarks, the parameter A2SC occurring 
in our analysis, which describes the effect of wd-pairing 
on the E/N of our strangelets, will not have exactly the 
same value as that describing 2SC pairing in bulk matter. 
Had we taken A2SC significantly smaller, our analysis of 
2SC pairing would not even be qualitatively reliable. 

We conclude that if the attraction between quarks is 
strong enough that the critical strangelets are in the 2SC 
phase (but not so strong as to favor the CFL phase for 
bulk neutral quark matter) then (Tcrit is increased by a 
factor which could be as large as « 2v^ if A2SC were 
tuned such that fi^^ -^ /x"^"^^"'' , but which is more 
typically smaller, of order 2. 



V. DISCUSSION 

If the surface tension in QCD of an interface between 
quark matter and vacuum is below a critical value Ccrit, 
and if the strange quark matter hypothesis holds, as- 
trophysicists may observe (or may be observing) strange 
quark stars that have crusts made of positively charged 
strangelets, with size of order the Debye length A^i, im- 
mersed in an electron gas. And, larger strangelets will 
be unstable to fission. In Sections |nl and IIIII we have 
developed a model-independent evaluation of cTcrit , along 
the way constructing profiles of strangelets taking into 
account both the energy benefit of charge separation and 
the energy costs introduced by the surface tension and 
Debye screening. Our result is given in H27|l . in terms of 
two parameters nq (the charge density of quark matter 
with /Xg = 0) and xq (the charge susceptibility of quark 
matter with ^g = 0) which are not currently known from 
first principles. Given any equation of state, whether 
from a model or from a future full QCD calculation, nq 
and XQ can be evaluated in terms of the quark number 
chemical potential /i and the strange quark mass ms. 

Our purpose in Section IIVI was to get a sense of the 
size of CTcrit, in MeV/fm^, using as a guide bag models 
for unpaired quark matter and quark matter in the 2SC 
color superconducting phase. A new twist in this section 
was the possibility that neutral bulk quark matter and 
stable strangelets could be in different phases (unpaired 
and 2SC respectively). We conclude from the model- 
dependent investigation in this section that it is easy to 
find parameters for which CTcrit ~ 1 — 3 MeV/fm^. And, 
by fine-tuning both nig and A2SC, it is possible to push 
CTcrit to as large as 5— 7 MeV/fm^. Although these conclu- 
sions are model-dependent, our results phrased in terms 
of ng and XQ as in Section llTTl can be applied in future to 
any model, for any phase, as our understanding of dense 
quark matter continues to improve. 

The surface tension ct is not known in QCD, but 
it has been calculated in the bag model for unpaired 



quark matter [SJI. It ranges from about 10 MeV/fm^ 
for ms = 100 - 120 MeV to about 4 MeV/fm^ for 
iTis = 250 MeV. If the interface has a thickness of or- 
der 1 fm, rather than being infinitely sharp as assumed 
in the bag model, the surface tension is almost certainly 
larger. It seems most likely that ct > CTcrit, making large 
strangelets stable against fission and giving strange quark 
stars, if they exist, fluid surfaces as in Ref. |a]. However, 
the combination of the bag model estimate for ct from 
Ref. [SJI and our model-dependent investigation in Sec- 
tion |W] leave open the possibility of the opposite conclu- 
sion. 

Further investigation of the properties of the strange 
star crusts that result from fragmentation and charge 
separation at strange star surfaces, as proposed in 
Ref. [3, is warranted. We have only analyzed the surface 
of the crust, where Poxt = and strangelets are spher- 
ical. If CT is less than CTcrit, we must then analyze the 
deeper layers of the crust. Preliminary numerical work 
indicates that Ag becomes less negative with increasing 
Poxt, but to date this ignores the alternative geometrical 
shapes of the strangelets that will occur below the sur- 
face of the crust. Furthermore, in the inner half of the 
crust, characterized by electron-filled "voids" embedded 
in quark matter rather than by strangelets embedded in 
an electron-gas, the curvature energy which we have ne- 
glected must be included. The curvature energy is pos- 
itive for strangelets, acting, like the surface tension, to 
suppress small strangelets. For voids, however, the curva- 
ture energy is negative, acting in opposition to the surface 
tension. If charge-separated crusts do occur on strange 
stars, it seems likely that they will prove more similar to 
conventional neutron star crusts than to the fluid surface 
previously considered for strange stars. It remains to be 
investigated, however, how similar to neutron star crusts 
they prove to be, and in particular whether the rich X- 
ray burst phenomenology observed in accreting compact 
stars can be consistent with a strange star crust. 

It remains unlikely that strange star crusts can give rise 
to pulsar glitches, since glitches require a crust within 
with charged nuggets (nuclei or strangelets would be 
equivalent) are immersed in both an electron gas and a 
superfluid. However, there is no superffuid in the strange 
star crust proposed in Ref. Q- However, if quark mat- 
ter in the crystalline color superconducting phase occurs 
somewhere within the core of a strange star, this could 
be the locus in which pulsar glitches originate |33 . l33| . 

Another way in which observations could settle some of 
the questions we have raised here would be the discovery 
of strangelets in cosmic rays, with a peak in the distri- 
bution of A and Z corresponding to the baryon num- 
ber Astabio of the most favored stable (or metastable) 
strangelet. If ct < CTcrit, and if strangelets occur in cos- 
mic rays, any strangelets whose initial baryon number 
(in whatever astrophysical collision produced them) was 
much larger than ^stable would fission into strangelets 
with A > ^stable- Strangelets with A only slightly bigger 
than Astabic will not fission because of the high energy en- 
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ergy per quark of the smallest strangelet produced. We 
therefore expect a peak in the strangelet distribution at 
A w Astabie , tailing off to no strangelets at all over some 
range of larger values. If the estimates of Ref. [33 of the 
flux of strangelets in cosmic rays incident on the earth 
are correct, the single-strangelet detection capability of 
the AMS-02 detector, scheduled to go into operation on 
the International Space Station in a few years, extends 
to strangelets with all ^ < lO"' [33. For two particular 
choices of model and model parameters in Section IIVI 
we found critical strangelets with A w 600 and Z w 70 
in one case and A « 680 and Z « 105 in the other 
case. More generally, we expect that if ct < CTcrit the sta- 
ble strangelets will have values of A lying between many 



hundred and a few thousand. If quark matter has such a 
small surface tension, then both strangelets and a peak 
in their size-distribution should be within the discovery 
reach of AMS-02. 
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